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Ab initio calculations within the density-functional theory formalism are performed to investigate 
the stability and electronic properties of fluorinated bilayer graphene (bilayer fluorographene). A 
comparison is made to previously investigated graphane, bilayer graphane and fluorographene. Bi- 
layer fluorographene is found to be a much more stable material than bilayer graphane. Its electronic 
band structure is similar to monolayer fluorographene but its electronic band gap is significantly 
larger (about 1 eV). We also calculate the effective masses around the Y- point for fluorographene 
and bilayer fluorographene and find that they are isotropic, in contrast to earlier reports. Further- 
more, it is found that bilayer fluorographene is almost as strong as graphene, as its 2D Young's 
modulus is approximately 300 iVra -1 . 

PACS numbers: 61.48. Gh, 68.43.-h, 68.43.Bc, 68.43. Fg, 73.21.Ac, 81.05.Uw 



I. INTRODUCTION 

Since the first reports on the successful isolation of sta- 
ble two-dimensional crystals consisting of a single atom 
layer by Novoselov et. al. in 2004, [TJ [2] researchers 
have been looking for ways to employ these new mate- 
rials in real-world applications. Of particular interest is 
graphene, a single layer derivative of graphite. Graphene 
attracts a lot of attention with its high crystal quality 
and very promising electronic properties. [3] 

Despite its unique characteristics which might be ex- 
ploited in future applications, (3HZ] there exists a seri- 
ous obstacle to use graphene in electronics as we know 
it today. Graphene is a zero-gap semiconductor and the 
absence of a band gap is protected by inversion and time- 
reversal symmetry. [5j Several approaches have been 
used or proposed to open and control a band gap in 
graphene and graphene derivatives: i) mechanical modifi- 
cation such as cutting the graphene sample into nanorib- 
bons, [HI [TO] ii) application of a potential difference over 
a bilayer of graphene [TTHT3] and, finally, iii) chemical 
modification which confines the p z electrons into covalent 
bonds. [T1HT5] The last approach includes the adsorption 
of radicals such as oxygen, hydrogen and fluorine atoms 
on the surface of graphene which form chemical bonds 
with the carbon atoms and change their hybridization 
from sp 2 to sp 3 . The chemical derivatives of graphene 
preserve the two-dimensional character of graphene but 
they have vastly different electronic properties. This 
makes them possibly more interesting materials for ap- 
plications in electronics than intrinsic graphene. 

Among the possible derivatives of graphene there 
are two that have attracted special attention, namely 
graphane and fluorographene (graphene fluoride), which 
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result from the hydrogenation and fluorination of 
graphene, respectively. They are theoretically predicted 
PH21)] and experimentally observed [HHH ITS ] 127 ] 125] 
to form crystalline materials in contrast to e.g. graphene 

oxide. [2EH3 

These new two-dimensional crystals are currently the 
subject of a considerable amount of experimental and 
theoretical work. But the process of chemical modifi- 
cation is not restricted to monolayer graphene. It has 
been proposed that this process can be expanded to bi- 
layer graphene as well. [31] Ab initio calculations showed 
that the weak van der Waals forces between the graphene 
layers are replaced by much stronger covalent bonds that 
stabilize the structure and that, at full coverage, a bilayer 
analogue of graphane is formed. The electronic struc- 
ture of intrinsic monolayer and bilayer graphane are very 
similar, |32j but their mechanical/elastic properties are 
expected to be substantially different. 

The difference in elastomechanical qualities can be an 
important issue. It was recently demonstrated |33j that 
the roughness of monolayer graphane surfaces is consid- 
erably larger than that of graphene. This increase of 
the size of the ripple formation can be explained by the 
lower stiffness of graphane [34] and its different vibra- 
tional properties as compared to graphene. |35j One way 
to reduce the increased roughness is to consider bilayers 
instead of single layers. As is well known from experi- 
ment the ripple formation in bilayer graphene is strongly 
reduced, i.e. they are removed by the intcrlaycr inter- 
action. [5S] A similar effect can be expected for bilayer 
graphane where the interlayer interaction is even more 
important. 

In this paper we investigate whether it is also possi- 
ble to fluorinate bilayer graphene. We perform ab ini- 
tio calculations to determine the stability and structural 
properties of a fluorinated bilayer of graphene, hereafter 
called bilayer fluorographene (for a better notion of the 
chemical composition and the structure see Fig. [I). We 
find that bilayer fluorographene is a much more stable 
compound than bilayer graphane, although more fluori- 
nation is needed to induce interlayer covalent C-C bonds. 
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The structural properties of bilayer fluorographene are in 
between those of monolayer fluorographene and diamond, 
other characteristics are closer to initial fluorographene. 

This paper is organized as follows: first we describe 
the computational details of our ab initio calculations, 
followed by an investigation of the stability and the 
formation conditions of interlayer covalent C-C bonds. 
Next, the overall stability and the geometrical proper- 
ties of fully fluorinated bilayer graphene (i.e. bilayer flu- 
orographene) are examined. These properties are sub- 
sequently compared with those of diamond and fluoro- 
graphene. To conclude, we investigate the electronic 
band structure and band gaps of single layer and bilayer 
fluorographene, together with the effective masses of the 
various possible charge carriers and elastic properties rep- 
resented by their 2D Young's moduli. 




FIG. 1: Top (left) and side (right) view of bilayer fluoro- 
graphene. The carbon atoms of the two layers are given by 
different colors (shades) for clarity and the fluorine atoms are 
in light green (open circles with the smallest diameter). 



III. RESULTS AND DISCUSSION 



A. Fluorination of bilayer graphene 



II. CALCULATIONS 



All our calculations were done within the density func- 
tional theory (DFT) formalism as implemented in the 
VASP package with usage of the local-density approx- 
imation (LDA) and the Perdew, Burke and Ernzerhof 
[37] generalized gradient approximation (GGA) for the 
exchange-correlation functional. We made use of the pro- 
jector augmented wave method |38| and a plane-wave ba- 
sis set with an energy cutoff of 500 cV. The relaxation of 
atomic positions was performed with forces smaller than 
0.01 eV A" 1 . 

Three types of supercells were used in our calculations: 
3x3 and 2x2 supercells to study the adsorption proper- 
ties of fluorine on a graphene bilayer for different concen- 
trations and configurations of fluorine, and a 1 x 1 unit 
cell for the calculation of the properties and electronic 
band structure of fully fluorinated graphene and bilayer 
graphene. 

The sampling of the Brillouin zone was done for the 
different supercells with the equivalent of a 24 x 24 x 1 
Monkhorst-Pack [39j fc-point grid for the monolayer or bi- 
layer graphene unit cell (containing 2 carbon atoms per 
layer). Spin polarization was not included in the cal- 
culations because fluorination is not expected to induce 
magnetism in graphene. [30] 

Since periodic boundary conditions were applied in all 
three dimensions the height of the supercell was set to 
20 A to include enough vacuum to minimize the inter- 
action between adjacent layers. Additionally we have 
performed a convergence test with respect to the planar 
size of the supercell. We obtain an accuracy in binding 
energies of less than 0.1 eV and in bond lenghts of less 
than 0.01 A with the chosen 3x3 supercell. 



The main objective of this paper is to investigate 
chemisorption of fluorine atoms on both sides of bilayer 
graphene. Hydrogenation and fluorination of the carbon 
atoms of graphene changes their hybridization from sp 2 
to sp 3 . In the case of hydrogenation of bilayer graphene, 
it was demonstrated that stable interlayer C-C bonds are 
induced at high concentrations of adsorbates.[2I] In this 
section, we perform a similar study for the case of fluo- 
rination and highlight the similarities and differences as 
compared to hydrogenated bilayer graphene. We make 
use of the local density approximation, because this ap- 
proximation gives a better description of the interlayer 
interaction in bilayer graphene over the generalized gra- 
dient approximation. [4"T] 




FIG. 2: Bilayer of AB stacked graphene with the 4 different 
sublattices indicated by different shades (colors). The A and 
B' sublattice are situated on top of each other as illustrated 
by the dotted lines. The carbon atoms at the corners of the 
light triangles indicate the 6 neighbors to be fluorinated. 

We only consider chemisorption on bilayer graphene 
with AB-stacking because this is the most energetically 
favorable way of stacking graphene layers. As indicated 
in Fig. [2j this gives rise to four sublattices (2 in each 
layer) of which two, namely the A and B' sublattice, are 
exactly aligned on top of each other while the sublattices 
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B and A', are not. Since the adsorbates under investiga- 
tion always chemisorb on top of a single C atom, there 
are two (inequivalent) adsorption sites, namely on the A 
(or B') sublattice and on the B (or A') sublattice. It 
was found in an earlier study of hydrogen adsorption on 
bilayer graphene that there is a difference in adsorption 
energy of approximately 0.03 eV in favor of the B sub- 
lattice. [3T] The H atoms are therefore likely to attach 
to those carbon atoms which have no direct neighbors in 
the other layer. For F atoms, we find the same preference 
for adsorption on the B sublattice. However, the differ- 
ence in adsorption energy for the two sublattices is only 
0.3 meV, which is about 100 times smaller than in the 
case of hydrogen. This is caused by the different nature 
of C-F bonds as compared to C-H bonds: C-H bonds 
are purely covalent, but C-F bonds have a partially ionic 
character because of the large difference in electronega- 
tivity between C and F. 

The formation of a chemical bond between an adsor- 
bate and a C atom of the bilayer of graphene induces 
a change in the hybridization of this C atom from sp 2 
to sp 3 . If the concentration of attached adsorbates on 
the preferred sublattice increases, the carbon atoms of 
the other sublattice also change their hybridization to 
sp 3 . This allows for the formation of covalent bonds be- 
tween the two graphene layers. These interlayer bonds 
are formed between carbon atoms of the A and B' sub- 
lattices that are pushed to one another because of the 
changed hybridization of their neighbors. 

For a single interlayer C-C bond there are six neighbors 
- three atoms in the B and three in the A' sublattice (Fig. 
[2|. However, it is not necessary for all these neighbors 
to change their hybridization in order for the interlayer 
bond to be formed, as in the case of hydrogen chemisorp- 
tion on bilayer graphene where only four of the neighbor- 
ing carbon atoms need to be hydrogenated to induce a 
covalent interlayer bond. |31j 

In order to study the formation of interlayer bonds 
upon fiuorination in more detail, we gradually increase 
the number of attached fluorine atoms on neighboring 
carbon atoms from 1 to 6. We completely relax the re- 
sulting systems and calculate the average and minimum 
distance between the C atoms of the two layers. The re- 
sults of these computations are summarized in Table [I] 
From this table, a gradual decrease of the interlayer dis- 
tance can be seen upon increasing the concentration of F 
atoms. At higher concentrations, however, a substantial 
jump in the minimum distance, d m i ni occurs when going 
from 4 to 5 fluorinated neighbors. This jump is about 
l A in size and clearly indicates the formation of an in- 
terlayer bond. It is thus necessary to fluorinate 5 of the 
6 neighbors to induce an interlayer bond, in contrast to 
hydrogenation where 4 hydrogenated neighbors suffice. 

Additional information about the stability of the inter- 
layer bond can be obtained by examining the formation 
and binding energy of the system at different concen- 
trations of adsorbed F atoms. We define the formation 
energy (Ef) as the energy per attached fluorine atom 



TABLE I: Interlayer chemical bond formation in a 3 x 3 
supercell. The binding energy per F atom (Et), the formation 
energy per F atom (Ef) and the average (d a v g ) and minimum 
distance (d m in) between the C atoms of the two layers are 
shown for different number of adsorbed fluorine atoms and 
different configurations. Presence of interlayer bond is marked 
as chemical bond. Energies are given in eV and distances in A. 



# F 


configuration 


E b 


Ej 


davg 




chem. 


atoms 












bond 


1 


B 


-2.525 


-0.819 


3.326 


3.256 


no 


2 


BB 


-2.405 


-0.700 


3.326 


3.190 


no 


2 


BA' 


-2.555 


-0.850 


3.213 


3.085 


no 


3 


BBB 


-2.290 


-0.585 


3.332 


3.126 


no 


3 


BBA' 


-2.487 


-0.782 


3.176 


2.979 


no 


4 


BBB A' 


-2.384 


-0.679 


3.150 


2.903 


no 


4 


BBA' A' 


-2.479 


-0.773 


3.061 


2.821 


no 


5 


BBB A' A' 


-2.450 


-0.745 


2.678 


1.757 


yes 


6 


BBB A' A' A' 


-2.560 


-0.854 2.598 


1.737 


yes 



with respect to intrinsic bilayer graphene and the di- 
atomic molecule F 2 . The binding energy (Ej,), on the 
other hand, is defined as the energy per fluorine atom 
(or CF pair) with respect to intrinsic bilayer graphene 
and atomic fluorine. Both energies are given in Table [T] 

The formation energy, Ef, is negative in all cases 
which means that all the configurations are stable against 
molecular desorption from the graphene surface. This 
should be contrasted to the same concentrations of hy- 
drogen atoms on the surface of bilayer graphene, where 
the energy becomes negative only in the case of almost 
fully hydrogenated bilayer graphene. [31| The value of 
the formation energies in the case of fiuorination is sig- 
nificantly lower (almost 1 eV) than in the case of hy- 
drogen chemisorption. This can be attributed to the 
weaker bond in the F 2 as compared to the H 2 molecule 
and is similar to the previous studied case of monolayer 
graphane and fluorographene.[23] The calculated bind- 
ing energies also show that it is energetically favorable 
for the fluorine atoms to attach themselves on both sides 
of the bilayer: for the same number of chemisorbed F 
atoms, the configuration in which these atoms are dis- 
tributed as equally as possible between the two layers 
is lower in energy and thus more stable. The fact that 
adsorbed F atoms on one side of the bilayer make it fa- 
vorable for other F atoms to adsorb on the other side, in- 
creases the chance of interlayer bond formation. Because 
these bonds are stable, we can imagine this process of fiu- 
orination and interlayer bond formation to continue until 
a fully covered bilayer of graphene fluoride is formed. 

To test if the aforementioned interlayer bond creation 
scenario is truly energetically favorable we considered a 
2x2 supercell for a gradual fiuorination with different 
configuration patterns. To be able to choose from the 
large number of different possible configurations we need 
some guidelines for further investigation. As a first guide- 
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line we distribute the F atoms equally on both sides of the 
bilayer because it is more energetically favorable as dis- 
cussed above. Four prevalent configurations are known 
for the fluorination of graphene: chair, boat, zigzag and 
armchair. The chair configuration can be readily ex- 
tended to the fluorination of bilayer graphene discussed 
above. The remaining three configurations share as a 
common feature that they contain dimers, as depicted in 
the inset of Fig. [3j Those dimers are used as the second 
guideline to reduce configuration space. The dimer a is 
composed of F atoms bonded with directly neighboring 
C atoms whereas the dimer b consists of F atoms bonded 
with distant C atoms on the opposite sides of the hexag- 
onal ring. For higher fluorine concentration (above two 
dimers, one at each side of the bilayer) dimers a and b are 
combined and create a trimer configuration, as depicted 
in the inset c of Fig. [3j 

The calculated results are given in Fig. [3j In this figure 
we present the formation energy per adatom for the dif- 
ferent fluorination pathways. These results show that for 
low concentrations of F atoms the dimer b configuration 
is more stable (carbon atoms from both sublattices are 
fluorinated and no interlayer bonds are created). How- 
ever, higher concentrations of F atoms are not feasible 
with the dimer configurations as the formation energy 
increases. But if only the B and A' sublattices are flu- 
orinated, the formation energy can decrease even more, 
finally creating bilayer fluorographene, with completely 
saturated C atoms as can be seen in Fig. [T] It is also 
important to point out that the total formation energy 
for all adsorbed F atoms (not shown in the figure) is a 
monotonically decreasing function of the F atom concen- 
tration and therefore the local minimum for dimer b will 
not result in an interruption of the adsorption. 

Our calculations also showed that for both dimer con- 
figurations no interlayer bonding was formed. The dimers 
a and b transform for higher coverage into the trimer 
configuration as indicated by the connection of the plots. 
For the adsorption on B and A' sublattices the interlayer 
bond is created already with only four F atoms. This ap- 
parently different result from the previous 3x3 supercell 
calculations is due to the choice of the size of supercell. 
Four adatoms in a 2 x 2 supercell correspond to a higher 
percentage of F per C atoms than five adatoms in the 
3x3 supercell. 

We also performed calculations for the case when the 
A and B' sublatices, lying on top of each other, were 
fluorinated. As already mentioned the preference for B 
(A') over A (B') sublattice for a single adatom is very 
small and unlikely to conduct the adsorption pattern. 
With increased F atom coverage the formation energy 
per adatom for fluorination of the A and B' sublattice 
was found to be an increasing function of the total cover- 
age. This results in a very unfavourable alternative over 
fluorination of B and A' sublattice. 

From the calculations on the F adsorption in a 2 x 2 su- 
percell we can conclude that the chair configuration is the 
lowest energy configuration, the only configuration that 



leads to interlayer bonding and in which a high fluorine 
concentration can be reached. 
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FIG. 3: Gradual fluorination of a 2 x 2 supercell of bilayer 
graphene. For increasing F concentration the formation en- 
ergy per fluorine atom is shown. The curves show different 
distributions of the F atoms - according to dimer configura- 
tion a (inset a) and b (inset b), only on sublattices B and A', 
and the placement of the F atoms only on sublattices A and 
B'. Dimers a and b can not be preserved for higher levels of 
fluorination, instead they transform into the trimer configu- 
ration (inset c). 



B. Properties of bilayer fluorographene 

In this section we examine the properties of bilayer 
fluorographene in more detail and compare them with 
those of monolayer fluorographene and diamond. For 
the study of the geometrical properties we include also 
GGA calculations because they give more accurate bond 
lengths and angles. A summary of the geometrical and 



electronic properties is given in Tables [TT] and III The 
results for monolayer fluorographene compare well with 
previous calculations on this system. [21] As can be ex- 
pected, the cell size, bond lengths and angles of bilayer 
fluorographene have values between those of monolayer 
fluorographene and diamond. Overall the geometrical 
properties are close to the one of bulk diamond due to 
the same hybridization of the carbon atoms. [25] 

The calculated C-C bond lengths are larger than the 
ones in diamond. This difference can be explained from 
a chemical point of view as a result of the depopulation 
of the bonding orbitals between the carbon atoms after 
substantial charge transfer to the F atom, which is the 
result of the difference in electronegativity between C and 
F atoms. 

Additionally we have observed the shortening of the C- 
F bond with increasing amount of fluorination on both 
sides of bilayer fluorographene. In a 3 x 3 supercell, the 
C-F bond length varies from 1.50 A for configuration B 
to 1.43 A for configuration BBBA'A'A'. As can be seen 
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from Tabic [TTJ the C-F bond length decreases even more 
to 1.38 A for fully ffuorinated bilayer graphene. The 
source of this behaviour can be attributed to the ionic 
character of the C-F bond and a gradual transformation 
from sp 2 to sp 3 hybridization of the C atoms. [4"0] 

TABLE II: Properties of single layer and bilayer fluoro- 
graphene: the unit cell length (a), the distances (d) and angles 
(8) between neighboring atoms. Distances are given in A and 
angles in deg. 





fluorographene 


bilayer 


diamond 








fluorographene 








LDA 


GGA 


LDA 


GGA 


LDA 


GGA 


a 


2.555 


2.596 


2.525 


2.563 


2.499 


2.527 


dc-c 


1.553 


1.576 


1.541 


1.563 


1.531 


1.547 


d C -c 


n/a 


n/a 


1.537 


1.554 


1.531 


1.547 


dc-F 


1.365 


1.382 


1.361 


1.377 


n/a 


n/a 


Occc 


110.7 


110.9 


110.0 


110.2 


109.5 


109.5 


8ccc 


n/a 


n/a 


108.9 


108.7 


109.5 


109.5 


Occf 


108.2 


108.0 


108.9 


108.7 


n/a 


n/a 



The electronic properties of monolayer and bilayer fluo- 
rographene are given in Table [TTl) we calculated the band 
gap of these materials together with the formation energy 
per atom (in contrast to previously used formation en- 
ergies per fluorine atom). For comparison, the values of 
these quantities are also given for diamond. 

The calculated band gap for monolayer fluorographene 
is in good agreement with previously published results 
when using GGA calculations. 42J Although the com- 
puted band gap value for fluorographene is close to ex- 
perimental measured values [H], this should be seen as 
a coincidence. LDA and GGA calculations are known to 
largely underestimate the value of the band gap. The cal- 
culated band gaps are in fact two times lower than more 
accurate results of many-electron GW calculations. [S3] 
This apparent discrepancy is attributed to the presence 
of a considerable amount of defects in the experimental 
samples which induce midgap states (similar to defected 
graphane 03]). 

The band gap of bilayer fluorographene is found to be 
larger than the one of monolayer fluorographene by ap- 
proximately 1 eV. This is different from the case of hy- 
drogenated graphene where monolayer graphane is found 
to have a slightly larger band gap than bilayer graphane. 
[3T] In fact, within LDA and GGA calculations, graphane 
has a larger band gap than monolayer fluorographene 
[3H US] but this observation does not apply to bilayer 
compounds where bilayer fluorographene surpasses bi- 
layer graphane in the size of the energy gap. Although 
our calculations are not accurate enough to provide the 
real band gap, it is probable that this difference in the 
size of the band gap between monolayer and bilayer fluo- 
rographene is qualitatively correct. This follows from the 
fact that LDA and GGA usually produce correct trends 
in the variation of the band gap among similar systems. 



So the band gap of bilayer fluorographene has been found 
to be intermediate between fluorographene and bulk di- 
amond. In this sense diamond can be considered as the 
limit of multilayer fluorographene. This tendency was 
observed before — for theoretically proposed graphite 
fluoride materials involving carbon atoms with sp 3 hy- 
bridization. |46j 

The absolute value of the formation energy of bilayer 
fluorographene is smaller than the one of the monolayer 
but still larger than the one of diamond. The main rea- 
son for the observed weakening of the formation energy 
is the drop of the ratio of the amount of fluorine atoms 
per carbon atom in going from monolayer to bilayer fluo- 
rographene. Overall, the qualitative image of the sta- 
bility of the fluorinated materials corresponds to that 
of graphane and bilayer graphane in terms of the afore- 
mentioned formation energy and the presence of covalent 
bonds between the graphene layers which stabilize the 
structure. Nevertheless there is a large quantitative dif- 
ference, the fluorinated materials are much more stable 
structures than the hydrogenated ones. [20] 

TABLE III: Electronic band gap E gap and formation energy 
per atom (Ef /atom) of monolayer and bilayer fluorographene 
and diamond using LDA and GGA calculations. All the en- 
ergies are given in eV. 





fluorographene 


bilayer 


diamond 






fluorographene 






LDA GGA 


LDA GGA 


LDA GGA 


Egap 


2.960 3.089 


3.951 4.040 


5.618 5.572 


Ef /atom 


-1.057 -0.944 


-0.722 -0.593 


-0.011 0.132 



The band structure of bilayer fluorographene is dis- 
played in Fig. [4] together with a diagram of the density 
of states. The depicted band structure is seen to be qual- 
itatively similar to the one of monolayer fluorographene 
[2"2"] with the size of the band gap as the only obvious 
difference. We also calculated the effective masses of the 
charge carriers around the T-point in order to find other 
differences or similarities between the two materials. 

The obtained effective masses of electrons, and heavy 
and light holes for both studied materials, can be found 
in Table IIVI It should be noted that these masses were 
calculated within the DFT formalism with usage of the 
LDA for the exchange-correlation functional which has 
been shown to give reasonable results. [4"7] 

The effective mass of charge carriers in crystalline ma- 
terials usually depends on the direction in reciprocal 
space. Therefore we have chosen two common directions 
for the hexagonal lattice at the T-point, namely: r — s- _K~ 
and r — ¥ M. The results for these two directions were 
found to be indistinguishable and so we can conclude 
that the effective masses at the T-point are isotropic. 
The direction independence for the effective masses that 
we observe contradicts previous results in which the ef- 
fective masses of monolayer fluorographene were found 
to be highly anisotropic. [35] Our statement about the 
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FIG. 4: The electronic band structure and density of states 
of bilayer fluorographene. The energies are relative to the 
Fermi level (Ef=0). 

isotropic character of the effective masses is supported 
by a direct plot of the energy value map around the T- 
point. In Fig. [5] we show a picture of the highest valence 
band of monolayer fluorographene. The isotropic char- 
acter of this band is clearly visible close to the T-point. 
Farther away from the T-point an anisotropy of the sur- 
face map is visible which can be attributed to fourth and 
higher order effects. However, these higher order effects 
do not induce anisotropy in the effective mass which is 
a second order effect. Similar results are found for the 
light hole and the electrons. 



ference in its values. We therefore conclude that both 
materials have similar electronic properties apart from a 
difference in the band gap of about 1 eV. 

TABLE IV: Effective masses of holes and electrons (in units 
of the electron mass) and the 2D Young's modulus, E, along 



i axes. E is expressed 


in Nm" 1 . 




fluorographene 


bilayer 




fluorographene 


m 


0.48 


0.50 


m hh 


1.13 


1.10 


mih 


0.41 


0.37 


E, r 


195 


284 


Ey 


197 


293 



In addition to the electronic properties we calculated 
the (2D) Young's modulus, E, for monolayer and bi- 
layer fluorographene which are also included in Table 
|IV| We followed the same approach as applied before in 
Ref. [22j The 2D Young's modulus of graphene is found 
to be 307 Nm^ 1 shifted from the experimental value, 
Eexp = 340 ± 50iVm _1 , 4§] and other theoretical values. 
[3"4"] The 2D Young's modulus of fluorographene is found 
to be 1/3 smaller in comparison to intrinsic graphene, 
whereas the obtained moduli for bilayer fluorographene 
reach almost the values of graphene, making these mate- 
rials very strong. The values unveil the isotropic charac- 
ter of the Young's modulus for both compounts. 



Energy levels (eV) 




FIG. 5: Energy levels of the highest valence band of mono- 
layer fluorographene in reciprocal space. The displayed region 
is a squared portition of the Brioullin zone with a side length 
0.22 of the T-M distance. The T-point is taken as the origin 
and the x-axis directed towards the M -point. Contour lines 
are displayed for a better notion of the symmetry. 



IV. CONCLUSIONS 

As an extension to previous work (22l |31] we studied 
the potentially interesting case of bilayer fluorographene. 
We demonstrated that fluorination of bilayer graphene 
results in more stable structures than hydrogenation. 
This can be clearly observed by a comparison of the for- 
mation energies of the final structures [5U] and is accen- 
tuated by the fact that the formation energy of partially 
fluorinated bilayer graphane is always negative in distinct 
contrast to partially hydrogenated bilayer graphene. The 
creation of interlayer chemical bonds occurs at higher 
amounts of fluorination as compared to hydrogenation. 

The calculated band gap of bilayer fluorographene 
shows a 30% increase over the one of bilayer graphane, 
and we also observed quantitatively significant differences 
between monolayer and bilayer fluorographene. From the 
value of the Young's modulus, we can conclude that bi- 
layer fluorographene is substantially stronger than mono- 
layer fluorographene and is almost as strong as graphene. 
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